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Abstract. We have derived elemental abundances of O, 
Na, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, Co, Ni as well as 
for a number of s-elements for 47 G and K dwarf, with 
[Me/H]> 0.1 dex. The selection of stars was based on 
their kinematics as well as on their uvby — /3-photometry. 
One sample of stars on rather eccentric orbits traces the 
chemical evolution interior to the solar orbit and another, 
on circular orbits, the evolution around the solar orbit. A 
few Extreme Population I stars were included in the latter 
sample. 

The stars have -0.1 dex < [Fe/H] < 0.42 dex. The 
spectroscopic [Fe/H] correlate well with the [Me/H] de- 
rived from uvby — (3 photometry. We find that the elemen- 
tal abundances of Mg, Al, Si, Ca, Ti, Cr and Ni all fol- 
low [Fe/H]. Our data put further constraints on models of 
galactic chemical evolution, in particular of Cr, Mn and Co 
which has not previously been studied for dwarf stars with 
[Me/H] > 0.1 dex. The increase in [Na/Fe] and [Al/Fe] as 
a function of [Fe/H] found previously by Edvardsson et 
al. (1993a) has been confirmed for [Na/Fe]. This upturn- 
ing relation, and the scatter around it, are shown not to be 
due to a mixture of populations with different mean dis- 
tances to the galactic centre. We do not confirm the same 
trend for aluminium, which is somewhat surprising since 
both these elements are thought to be produced in the 
same environments in the pre-supernova stars. Nor have 
we been able to trace any tendency for relative abundances 
of O, Si, and Ti relative to Fe to vary with the stellar ve- 
locities, i.e. the stars present mean distance to the galactic 
centre. These results imply that there is no significant dif- 
ference in the chemical evolution of the different stellar 
populations for stars with [Me/H]> 0.1 dex. We find that 
[O/Fe] continue to decline with increasing [Fe/H] and that 
oxygen and europium correlate well. However [Si/Fe] and 
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[Ca/Fe] seem to stay constant. A real ("cosmic") scatter in 
[Ti/Fe] at given [Fe/H] is suggested as well as a decreas- 
ing abundance of the s-elements relative to iron for the 
most metal-rich dwarf stars. We discuss our results in the 
context of recent models of galactic chemical evolution. 

In our sample we have included a few very metal rich 
stars, sometimes called SMR (super metal rich) stars. We 
find these stars to be among the most iron-rich in our 
sample but far from as metal-rich as indicated by their 
photometric metallicities. SMR stars on highly eccentric 
orbits, alleged to trace the evolution of the chemical evo- 
lution in the galactic Bulge, have previously been found 
overabundant in O, Mg and Si. We have included three 
such stars from the study by Barbuy & Grenon (1990). 
We find them to be less metal rich and the other elemen- 
tal abundances remain puzzling. 

Detailed spectroscopic abundance analyses of K dwarf 
stars are rare. Our study includes 5 K dwarf stars and has 
revealed what appears to be a striking example of overi- 
onization. The overionization is especially prominent for 
Ca, Cr and Fe. The origin of this apparent overionization 
is not clear and we discuss different explanations in some 
detail. 

Key words: Galaxy: abundances - evolution - Stars: 
abundances - late-type - fundamental parameters 



1. Introduction 

The time that F and G dwarf stars spend on the main- 
sequence span a range from 10 9 to several times 10 10 years. 
This means that such stars may be used as tracers of the 
chemical and dynamical evolution of the Galaxy; in fact, 
a combination of chemical and kinematical data is a very 
powerful tool for studying the galactic chemical evolution, 
cf. e.g. Edvardsson et al. (1993a) and Wyse & Gilmore 
(1995). 

Spectroscopic abundance analysis of stars have now be- 
come accurate enough to admit determinations of rather 
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small 0.10 dex) relative abundance differences in differ- 
ential studies of stellar samples of large size. These devel- 
opments make it possible to explore the galactic chemical 
evolution in considerable detail. 

Edvardsson et al. (1993a) analysed 189 F and G dwarf 
stars, with —1.1 dex < [Fe/H] < 0.25 dex. Accurate abun- 
dances were determined for a number of key elements, O, 
Na, Mg, Al, Si, Ca, Ti, Fe and Ni, as well as a number 
of s-process elements. The abundance results combined 
with accurate velocity data enabled a detailed study of 
the chemical evolution of kinematically distinct popula- 
tions. Extensions of this study were made by Tomkin et 
al. (1995) and Woolf et al. (1995) who measured carbon 
and europium abundances, respectively, for about half of 
the stars in the Edvardsson et al. (1993a) sample. 

The study by Edvardsson et al. (1993a) raised a num- 
ber of new questions concerning the most metal-rich stars 
in the galactic disk, not the least concerning the build- 
up of sodium, magnesium and aluminium. Sodium, alu- 
minium and, possibly, magnesium relative iron vs [Fe/H] 
showed an increase for [Fe/H] > 0.0 dex (cf. Edvardsson et 
al. 1993a Fig. 15a-l). Are these "upturns" real? Large star- 
to-star scatter was also encountered for the abundance of 
certain elements: magnesium, aluminium and titanium, 
relative to iron at a given [Fe/H]. Could this scatter be 
reduced by using more/better abundance criteria or is the 
scatter intrinsic to the stellar population? One suggestion 
was that "upturns" and scatter could be due to a mix- 
ing of populations with different ages and with different 
mean distances from the galactic centre, e.g. a mixture 
of old metal-rich stars, more concentrated to the centre 
and young Extreme Population I stars on solar like orbits. 
We have therefore studied a sample of 47 metal-rich dwarf 
stars, with photometric metallicities [Me/H] > 0.10 dex, 
chosen to represent different mean perigalactic distances 
and presumably different ages. 

The paper is organized as follows: in Sect. 2 and 3, 
we describe the selection criteria of the stellar sample, the 
observations and reductions, Sect. 4 contains a descrip- 
tion of the analysis while the errors are discussed in detail 
in Sect. 5. Abundance results and their interpretation in 
terms of models of galactic chemical evolution are pre- 
sented in Sect. 6 and, finally, Sect. 7 contains a summary 
and discussion. 

2. Selection of stars 

The stars have been selected to include both stars from 
the old metal-rich disk population, as well as more recently 
formed Extreme Population I dwarf stars. The sample of 
stars was confined photometrically by 0™38 < b—y < 0™63 
and with 8c\ and Ami such that with the calibrations of 
Edvardsson et al. (1993a) and Olsen (1984), solar-type 
dwarfs and subgiants are singled out with 4500K < T c g < 
6000 K, 4.0 < logg < 4.6 and [Me/H] > 0.10 dex. Magni- 
tude and metallicity data were taken from the catalogues 



by Olsen (1983, 1993, 1994 and priv. comm.). These were 
combined with accurate velocity data (UVW) calculated 
by Olsen from available proper motions and unpublished 
CORAVEL radial velocities by Andersen, Mayor & Nord- 
strom, in the selection of programme stars. 

From the catalogue we seleted 42 G and K dwarf stars 
with [Me/H] > 0.10 for our study. 13 of these stars have 
Vlsr < -50 km/s and/or Q LS r > 60 km/s, where Vlsr 
is the stellar space velocity component in the direction of 
the rotation of the Galaxy relative to the local standard of 
rest (LSR) and Qlsr is the total space velocity relative to 
LSR. These stars were selected to represent the chemical 
evolution galactic regions different from those represented 
by the solar orbit, in particular located more closely to the 
galactic centre. As shown in Edvardsson et al. (1993b), 
Fig. 1, Vlsr < —50 km/s singles out stars with mean 
perigalactica less than 7 kpc. The 17 stars with Qlsr < 30 
km/s represent the chemical evolution in the gas close to 
the solar orbit. 

Since we also wished to explore the oxygen abundances 
in metal-rich disk dwarf stars on highly eccentric orbits we 
added 3 additional stars from the sample of stars for which 
Barbuy & Grenon (1990) found abnormally high [O/Fe] 
ratios (HD37986, HD77338, HD87007). Moreover, we in- 
cluded 2 stars from Schuster & Nisscn (1988) (HD175518 
and HD182572). 

The total final sample contains 47 stars and has 
3 in common with Edvardsson et al. (1993a): HD30562, 
HD67228 and HD144585. 

3. Observations and reductions 

3.1. Observations 

The observations were carried out with the 2.7 m telescope 
of the McDonald Observatory, University of Texas, during 
two observing runs, March 1 — 9, 1994 and April 28 - 
May 2, 1994, using the 2d-coude cross dispersion echcllc 
spectrometer (Tull et al. 1995). 

At the time of the observations the new 2kx2k chip, 
which covers most of the visual spectrum of the observed 
star in one exposure, was not yet available at the spec- 
trometer. The CCD-chip used has 800x864 pixels and 
covers approximately 18 echelle orders of length 25—35 
A. Settings of the detector for the observations were se- 
lected to enable reliable analysis of the spectral lines 
from the elements oxygen, sodium, magnesium and alu- 
minium, silicon, calcium, titanium, as well as chromium, 
iron, nickel, yttrium and europium. Many sources were 
utilized to make the selection of lines optimized for the 
current project, the most important being Edvardsson et 
al. (1993a), Morell (1994), Mackle et al. (1975), Ruland 
et al. (1981). The Arcturus Atlas by Griffin (1968) and 
the Solar Spectrum Atlas by Delbouille et al. (1973), were 
carefully inspected, in order to avoid blends and define 
suitable continuum regions. Altogether 3 different settings 
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of the CCD in the focal plane of the spectrometer camera 
were used. 

The nominal spectral resolution (slit width 0.25 /im) 
of ~ 100 000, was verified by observations of a thorium- 
argon lamp and by measuring telluric lines. This relatively 
high resolution was judged to be important in view of 
the relatively crowded spectra, in particular for the cooler 
stars. 

3.2. Reductions of stellar spectra 

The reductions were carried out with standard IRAFQ 
packages for reduction and extraction of cchelle spectra. 
The procedure adopted was: all raw frames were first cor- 
rected for read-out noise and bias, the flat fields added to- 
gether, all object frames were trimmed, the new flat field 
frames were normalized, the stellar frames were divided by 
the normalized flat field, corrected for scattered light and 
the one-dimensional stellar spectra extracted. Finally, the 
spectrum of the comparison lamp, a thorium-argon lamp, 
was reduced, extracted and the lines in it identified. The 
comparison spectra were used to carry out the wavelength 
calibration of the stellar spectra. A Legendre polynomial 
was fitted to each stellar spectrum to define the contin- 
uum, using the CONTINUUM task in IRAF. 

The results of this procedure for defining the contin- 
uum were inspected visually and judged to give good con- 
tinua for all the lines used in the abundance analysis; how- 
ever, for the strong Cal line at 6162 A, used for checking 
surface gravities, this definition of the continuum was not 
satisfactory for all stars. Instead, a blaze function was con- 
structed for each night by fitting a Legendre polynomial to 
the spectrum of a B star observed during the night. The 
stellar spectra were then divided by this blaze function. 
Then a low order Legendre polynomial was fitted to the 
stellar spectrum to rectify it. The change in continuum 
was negligible for half of the stars, as compared with the 
result of using CONTINUUM; this also includes half of 
the K dwarf stars. 

Our spectra are generally of high quality, reaching S /N 
of ~ 200 in most cases. Some problems with fringes in 
the red were encountered. The measurements of spectral 
line equivalent widths were done using the SPLOT task in 
IRAF. Karin Eriksson is gratefully thanked for carrying 
out almost half of the measurements. Care was exercised 
to guarantee that the measurements agreed well between 
the two operators. 



1 IRAF is distributed by National Optical Astronomy Ob- 
servatories, operated by the Association of Universities for Re- 
search in Astronomy, Inc., under contract with the National 
Science Foundation, USA. 



3.3. Solar observations and reductions 

The solar observations were performed by illuminating the 
slit of the spectrograph by sky light each afternoon. The 
solar spectra have ~ 2 times higher S/N ratios than the 
stellar spectra. The reductions and measurements of the 
solar spectra were made in the same way as for the stellar 
spectra. 

The solar observations were used to determine "astro- 
physical" log gf- values. For lines, denoted by K in Table 
||, shifted by the difference in radial velocities between the 
stars and the Sun outside our solar spectral recordings, 
but judged to be of interest for the analyses, solar equiva- 
lent widths were instead obtained from the solar spectrum 
atlas by Kurucz et al. (1984). 

4. Analysis 

We have performed a standard Local Thermodynamic 
Equilibrium (LTE) analysis, strictly differential with re- 
spect to the Sun, to derive chemical abundances from the 
measured equivalent widths. 

4-1. Model atmospheres 

We used the MARCS program, first described by Gustafs- 
son et al. (1975), to generate the model atmospheres. Since 
then, the program has been further developed in various 
ways and updated in order to handle the line blanketing 
of millions of absorption lines more accurately, Asplund 
et al. (1997). The following assumptions enter into the 
calculation of the models: the atmosphere is assumed to 
be plane-parallel and in hydrostatic equilibrium, the total 
flux (including mixing- length convection) is constant, the 
source function is described by the Planck function at the 
local temperature with a scattering term, the populations 
of different excitation levels and ionization stages are gov- 
erned by LTE. Since the analysis is differential relative 
to the Sun we have also used a solar model atmosphere 
calculated with the same program as the stellar models, 
in spite of the fact that the empirically derived Holweger- 
Miiller model better reproduces the solar observed limb 
darkening; see Blackwell et al. (1995) for a discussion of 
this. 

4-2. Fundamental parameters of model atmospheres 

The effective temperature and surface gravity for each star 
were derived from uvby — (3 photometry from Olsen (1983, 
1993, 1994 and priv. comm.). We have used an extension 
of the calibration in Edvardsson et al. (1993a) and the cal- 
ibration in Olsen (1984). The calibration by Edvardsson 
et al. (1993a) is valid for the parameter space 2.58 < (3 < 
2.72 and 0.04 < (da + 0.5 ■ 6mi) < 0.16. This translates 
roughly to 5600 K < T cff < 7000 K, 3.8 < logg < 4.5 and 
metallicity -1 dex < [Me/H] < +0.3 dex. The Edvards- 
son et al. (1993a) calibration was preferred for all stars to 
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Table 1. Stellar parameters of the observed stars. The fifth column gives the spectral classification according to 
01scn(1983, 1993, 1994 and priv. comm.), except where indicated 1) from Eggcn (1960), and 2) spectral class from 
Hofflcit & Jascheck (1982); the sixth and seventh columns give effective temperature and surface gravity as derived 
from photometry; the eight and ninth give surface gravities as derived from spectra, wings of Cai, and parallaxes, 
respectively; the tenth and eleventh give iron abundances as determined from our spectral analysis and mctallicities as 
determined from photometry; the following four columns contain the space velocity data for the stars, in km/s relative 
to the LSR, from Olsen (priv. comm.). In the last column the stellar space velocity is given. 
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3 All stellar parameters from Barbuy & Grenon (1990) 

4 Barbuy, private communication 
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Fig. 1. Surface gravities derived from parallaxes vs sur- 
face gravities derived from photometry. The one-to-one 
relation is indicated by the dotted line. The error bars re- 
flect the errors in the parallaxes. Note the different scales 
on x- and y-axis. 



which it is applicable. For the three stars previously stud- 
ied by Barbuy & Grenon (1990) uvby — (3 photometry is 
not available, and we used the stellar parameters cited by 
them. 



4-3. Surface gravity 

Surface gravities were derived from the c\ index, which 
primarily measures the Balmer jump and is sensitive to 
the surface gravity in solar type dwarf stars. In addition 
to surface gravities estimated from photometry we have 
also used the wings of the strong calcium line at 6162 A 
(see Blackwell & Willis (1977) and Edvardsson (1988)), to 
derive surface gravities. The results are given in Table § 
For these results we are indebted to Matthias Palmer and 
Mikael Nilsson, who carried out these laborious determi- 
nations. Edvardsson et al. (1993a) estimated the error in 
surface gravity determined from their calibration of the c\ 
index to 0.2 dex. Olsen (1984) quotes a similar accuracy 
for his calibration. For the stars with trigonometric par- 
allaxes in van Altena et al. (1991) we have also estimated 
surface gravities from parallaxes. This exercise was meant 
to be a consistency check on the photometric and spectro- 
scopic determinations of the surface gravities. Assuming 
that the stars have masses of 0.8 M© and interpolating the 
bolometric corrections in the table given in Allen (1973) 
we derived, from the parallaxes, the surface gravities (see 
e.g. Gustafsson et al. 1974) shown in Fig. [I] and given in 
Table |l]0. 



2 With the release of the Hipparcos parallaxes this can now 
be addressed in greater detail and for most of our stars. The 
results will be published elsewhere. 



The agreement between surface gravities determined 
with the different methods is in general good. We note, 
however, that there is in the mean an offset in the surface 
gravities determined from the wings of the strong Ca I line 
as compared to the surface gravities determined from pho- 
tometry, see Table [l], of about +0.27 dex. The reason for 
this offset is not clear but, a corresponding uncertainty in 
surface gravity is of minor significance for the abundance 
results, see Table |[ A tendency for the trigonometric log g 
determiantions to agree closer with the photometric than 
with the spectroscopic ones may also be traced. It is in- 
teresting to note that the K dwarf stars in our study have 
logg derived from parallaxes that are in good agreement 
with those determined from photometry. 

On the basis of these comparisons we found no reason 
to change the surface gravities to be used in the abundance 
analysis but kept those determined from photometry. 

4-. 4- Microturbulence parameters 

We have determined microturbulence parameters, £t, from 
both the Ca I lines and the Fe I lines for 12 of the stars with 
enough calcium and iron lines measured. For the determi- 
nation from the Cai lines, abundances from the individual 
lines were derived with £ t ranging from 0.10 to 1.90 km/s. 
The microturbulence parameter for each star was then 
determined as the £t value which gave the smallest abun- 
dance scatter (the inflexion point, see e.g. Smith 1981). 
The microturbulence parameters were also determined by 
plotting the iron abundance versus the reduced equivalent 
widths, \og(W\/X), derived for each line. If the correct mi- 
croturbulence parameter is used the slope of a fit to the 
data points should be zero. Both methods agreed well for 
these 12 stars. 

Using the values obtained we derived a relation be- 
tween effective temperature, surface gravity and the mi- 
croturbulence parameter, which was then used to deter- 
mine microturbulence parameters for the rest of the stars: 
£ t = 4.5 • 10~ 4 T c ff — 0.31 • log g and predicts, with T e g and 
logg adopted, £ t values to an accuracy of ±0.1 km/s for 
the 12 stars. This relation is valid for 3.93 < \ogg < 4.50 
and 5530 K < T e ft < 6036 K. For a few stars just out- 
side the validity range we extrapolated the relation to cal- 
culate approximative microturbulence parameters. This 
procedure was checked to yield consistent abundances 
for individual iron lines of different strengths. For five 
stars a somewhat lower microturbulence parameter was 
preferred (HD36130, HD134474, HD171999, HD182572, 
HD187055). Our relation for microturbulence parameters 
yield slightly higher values than the relation presented by 
Edvardsson et al. (1993a). 

For the K dwarf stars none of the described methods 
seemed to yield definite values for the microturbulence 
parameter. A microturbulence parameter of 1.0 km/s was 
adopted for these stars. 
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Table 2. Column 1: Wavelength as quoted in Moore et al. (1966). Column 2 : Excitation energy of the lower level in- 
volved in the transition. Column 3: Astrophysically derived oscillator strengths (based on the solar equivalent widths). 
Column 4: Correction factor to Unsold's damping constant. Column 5: Radiative damping constant. Column 6: Equiv- 
alent width as measured in the Sun. If nothing else is indicated in the 7th column observations of the daylight sky 
with the same spectrograph were used. Column 7: Lines used to check the excitation equilibrium are denoted by t and 
lines for which the solar equivalent width was measured from the Kurucz et al. (1984) Solar Flux Atlas by a K, see 
Sect. B.3. Cm and Fen lines excluded from the analysis of the K dwarfs stars due to possible blends are denoted by 
+ , see Sect. 5.6. For each element the solar logarithmic abundance is given according to in Anders & Grevesse (1989), 
except for iron for which the value 7.51 was used, see Holweger et al. (1991) and Biemont et al. (1991). 
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4-5. Atomic line data 

The oscillator strengths were determined by requiring 
that the abundances calculated for the solar model 
(T cff =5780K, log 5 = 4.44, [Me/H]=0.00, 6=1-00) should 
reproduce the observed equivalent widths of the solar spec- 
trum. The resulting log gf values are given in Table § 

Different line broadening mechanisms, van der Waals 
damping, radiation damping, thermal Dopplcr broadening 
and microturbulence were considered in the calculations of 
equivalent widths and abundances. Enhancement factors 
for the van der Waals damping were compiled from the lit- 
erature. For iron lines values from Hannaford et al. (1992) 
and Holweger et al. (1991) were used, for calcium values 
from Smith (1981) and references therein, for sodium val- 



ues from Holweger (1971) and for scandium values from 
Neuforge (1992). For the remaining lines a correction fac- 
tor of 2.5 was adopted to the classical Unsold value, ac- 
cording to Mackle et al. (1975). For the (unimportant) 
radiation damping parameter values from Kurucz (1989) 
were adopted for lines from calcium through nickel. 

5. Errors in resulting abundances 

5.1. Measurement of equivalent widths and fitting of 
continua 

If we assume that the true continuum level is not wildly 
different from the fitted continuum the error in measured 
equivalent width can at the most be as large, in difficult 
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Table 3. Effects on abundance estimates for a number of elements from changes in fundamental parameters of the 
model atmosphere. The resulting abundances are shown for 3 stars: HD72946 (T o ff/log#/[Fe/H])=(5911/4.4/0.24), 
HD103932 (4510/4.58/0.16), and HD110010 (5965/4.08/0.35). The first line for each star gives [X/H], X being the 
element indicated in the header, derived using a model atmosphere constructed with the stellar parameters adopted in 
this study. The following lines contain the differences when effective temperature or surface gravity has been changed 
as indicated in the first column. Note that the changes in parameters are different, and larger, for HD103932. 



ID / A 


[Ol] 


Nai 


Mgi 


All 


Si 


Cai 


Til 


Cr I 


Cm 


Fei 


Fen 


Nil 


Eun 


HD72946 

+100 K 
+0.2 dex 


0.16 
+0.01 
+0.09 


0.22 
+0.04 
-0.03 


0.24 
+0.04 
-0.01 


0.23 
+0.06 
-0.01 


0.13 
+0.01 
-0.01 


0.22 
+0.06 
-0.02 


0.38 
+0.08 
0.00 


0.24 
+0.06 
0.00 


0.09 
-0.03 
+0.08 


0.24 
+0.05 
-0.01 


0.00 
+0.01 
+0.12 


0.18 
+0.05 
+0.01 


0.09 
0.00 
+0.08 


HD103932 

+200 K 
+0.4 dex 


0.28 
+0.05 
+0.16 




0.30 
-0.08 
+0.02 


0.14 
+0.12 
-0.04 


0.27 
-0.18 
+0.10 


-0.12 
+0.18 
-0.04 


+0.03 
+0.24 
-0.02 


0.01 
+0.12 
+0.02 


0.44 
-0.18 
+0.19 


0.14 
-0.04 
+0.07 


0.57 
-0.32 
+0.21 


0.41 
-0.05 
+0.09 


0.32 
-0.01 
+0.15 


HD110010 

+100 K 
+0.2 dex 


-0.04 
+0.01 
+0.09 


0.53 
+0.05 
-0.02 


0.48 
+0.04 
-0.01 


0.40 
+0.04 
-0.02 


0.51 
+0.01 
-0.01 


0.35 
+0.06 
-0.01 


0.29 
+0.08 
0.00 


0.30 
+0.06 
0.00 


0.22 
+0.03 
+0.08 


0.35 
+0.06 
-0.01 


0.33 
-0.05 
-0.03 


0.31 
+0.05 
0.00 


0.16 
0.00 
+0.09 



cases, as 2 mA for weak lines, i.e.~ 20%, and 4-5 mA 
for stronger lines, i.e.~ 10%. This translates to typically 
0.08 dex in the resulting abundance derived from a line 
of strength 10 mA and 0.04 dex for a line of strength 50 
mA. Apart from the effects of the continuum errors and 
blends, the error in derived abundances due to errors in the 
measurement of the equivalent width of a line is negligible 
(< 0.01 dex). In general, the lines and continua are, due 
to the high S/N and high resolution, easy to fit and the 
errors given above should be regarded as upper limits. 



5.2. Oscillator strengths 

The oscillator strengths, derived from the observed solar 
spectrum, can be affected by misidentification, by unde- 
tected blends and by errors in continuum fitting and mea- 
surements of equivalent widths in the solar spectrum. As 
in the stellar spectra, location of the continuum is a much 
larger source of error than the actual measurement of a 
line. (Note, however, that the solar spectra have higher 
S/N, usually ~ 400, and thus, identification of the contin- 
uum becomes easier as well as identification of lines. The 
Sun is also more metal-poor than the programme stars, 
which makes identification of the continuum in the Sun 
easier.) Using the results in Sect. 5.1 we find that errors 



in logg/-values may be as large as 0.08 dex, but a more 
characteristic number is 0.04 dex. 

Since we perform a purely differential analysis errors 
due to misplaced continua, neglected blends, etc. should 
partly cancel in the first approximation and not affect the 
resulting differential abundances very much, as long as we 
study stars similar to the Sun. 



5.3. Blends 

In the selection of lines we have carefully avoided all lines 
that could be subject to blending with nearby lines as 
given in Moore et al. (1966). For ions with several lines 
measured we have also looked for lines which produces ab- 
normally high abundances as compared with the majority 
of the lines. This led us to exclude three Fell lines from 
our final analyses: 6383.71, 6383.45, and 6627.32 A. 

5.4. Fundamental parameters of the model atmospheres 

Edvardsson et al. (1993a) estimate the error in the effec- 
tive temperature due to errors in b — y to range from -50 
K to +100 K and the corresponding error in logg to be 
±0.2 dex. 

The effects of such errors in T c g and logg are ex- 
emplified in Table g|. As expected, abundances derived 
from ions are most sensitive to changes in surface gravity 
while abundances derived from atoms are most sensitive 
to changes in effective temperature. In general, errors in 
derived abundances are smaller than 0.10 dex for atoms 
when varying the effective temperature by ±100 K and 
less than 0.02 dex when the surface gravity is varied by 
±0.20 dex; they are smaller than 0.02 dex for ions when 
the effective temperature is varied and less than 0.10 dex 
when the surface gravity is varied. 

5.5. Photometric versus spectroscopic metallicities 

For dwarf stars that are significantly more metal- rich than 
the Sun (e.g. [Me/H] sa 0.2 dex), the metallicity used in the 
calculation of the model atmosphere is important, since 
it governs the line blanketing and thus the temperature 
structure of the model atmosphere. We may expect that 
final derived abundances are sensitive to this parameter. 
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Fig. 2. Metallicitics determined from photometry vs iron 
abundances derived in our spectral analysis. The one-to- 
one relation is indicated by a dotted line and a least-square 
fit to the data points, taking the error in [Fe/H] to be 
±0.02 and the error in the photometric metallicities to 
be ±0.f, is also shown, dashed line. The x 2 -proabability 
for this fit is 0.92. A fit made without taking the errors 
into account is shown by a long-dashed line. The stars 
from Barbuy & Grenon (1990), •, and HD171999A and 
HD175518, denoted by o symbols, were excluded from the 



fit (see Sect. 5.5) 



Following Edvardsson et al. (1993a) we decreased (and 
in a few cases increased) the metallicities to the values 
derived for [Fe/H] in a preliminary abundance analysis 
and reiterated the abundance determination. (We have 
determined iron abundances from lines arising from both 
Fel and Fen. For Fei usually more than 30 lines were 
analysed and for Fe n three to four lines. The formal error 
in the mean [Fe/H] derived from Fei for a certain star 
is usually smaller than 0.02 dex.) In the mean we find 
that we had to reduce the abundances by 0.07 dex from 
the initial photometric values, with a spread of 0.01 dex, 
Table |l|. Apparently, our photometric metallicities tend to 
overestimate the metal content in these metal-rich stars 
as compared to the iron abundance derived from spectral 
abundance analysis, Fig.|2|. 

Usually, the differences between the metallicities es- 
timated from photometry, are close enough to those de- 
rived from spectroscopy that we have not found it nec- 
essary to reiterate the determination of effective temper- 
ature and log g. (The determination of metallicity was, 
however, changed according to the spectroscopic result so 
that the final model used in the analysis had [Me/H] con- 
sistent with the resulting spectroscopic [Fe/H].) However, 
for some of the stars the photometry indicates rather ex- 
treme metallicities compared with the spectroscopic iron 
abundances. For HD171999A we have only measured 6 Fe I 
equivalent widths (since th is s tar was only observed with 
one CCD setting, see Sect. iA) and thus the spectroscop- 



Fig. 3. Copper abundances derived from the lines at 5220 
A, • symbols, and 7933 A, □ symbols. The widths of these 
lines are 16 and 36 mA, respectively, in the solar flux spec- 
trum. 



ically determined iron abundance is not as good as for the 
other stars. However, we note that the line-to-line scatter 
is small, 0.03 dex. For HD175518 it is questionable if its 
photometric metallicity is realistic. In any large catalogue 
there will always be a few stars with 2 — 3cr errors in the 
photometry. Since spectroscopic iron abundances indicate 
a much lower metallicity this is most probably an example 
of that. 

We have studied, for HD175518, the effects on derived 
abundances if [Me/H] is lowered to 0.2 dex as indicated 
by the spectroscopy, thus affecting the estimates of the 
rest of the fundamental parameters. The lower metallicity 
implies a lower effective temperature; [Me/H] = 0.22 dex 
corresponds to a decrease of T c g by ~200 K. From Table 
U we see that most elements will change by ~ 0.1 dex 
and thus the star will mainly move horizontally, by ~ 0.1 
dex, in our relative abundance diagrams. HD175518 is an 
extreme case in our sample and the abundances of the 
other stars and general abundance trends for the whole 
sample should not be affected by comparable amounts. 

5. 6. Effects of hyperfine structure 

Some of the elements analysed are well known to be sub- 
ject to hyperfine structure. We have not taken this into 
account when deriving the chemical abundances. Instead, 
for those elements in particular we have aimed at selecting 
weak enough lines, so that the neglect of hyperfine struc- 
ture in the abundance analysis should not affect the cal- 
culated equivalent widths and thus not the derived abun- 
dances. 



Copper. Figure y shows copper abundances derived from 
the two lines used in our study as functions of equiva- 
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Fig. 4. Two examples of how the slope coefBcients in Fig. 
H were obtained. Examples are for a. HD91204 and b. 
HD 125968. The dashed lines represent linear least-square 
linear fits. 

lent width. Our data clearly show that the stronger cop- 
per line is subject to hyperfine structure and that these 
lines should be analysed using synthetic spectroscopy tak- 
ing the hyperfine structure splitting into account. Results 
by Summers (1994) may suggest departures from LTE in 
the population of levels in the copper atom. We do not, 
however, have enough data to make further empirical in- 
vestigations of such departures from LTE. We omit copper 
abundances from the following discussion. 

Manganese and Cobalt. The manganese and cobalt lines 
used in this study are not saturated. When plotting abun- 
dances derived from each line as a function of equivalent 
width no distinct pattern was found; indicating that the 
omission of hyperfine structure in the treatment of the 
lines is not problematic. In many stars lines with equiv- 
alent widths of about 20 and 60 mA, respectively, gave 
manganese abundances that are in excellent mutual agree- 
ment. 

r- and s-process elements. The lines used to derive abun- 
dances for the heavy s- and r-process elements are suffi- 
ciently weak in our programme stars to be safely used as 
abundance criteria in spite of being subject to hyperfine 
structure. Among the light s-process elements we note, 
however, that only Yn, and maybe Yi, have lines strong 
enough and secure enough that derived abundances can 
be used with confidence. The log ^/-values for zirconium 
are very uncertain because of the faintness of the lines in 
the Sun. 

5.7. Checks on Boltzmann and Saha equilibria 

Deviations from the Boltzmann excitation equilibrium, 
which might suggest an error in the effective tempera- 
ture, can be traced by studying abundances derived from 
individual lines as a function of the lower excitation po- 
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Fig. 5. The slope coefficients from the excitation energy 
balance diagram for each star plotted vs effective temper- 
ature, surface gravity and spectroscopically derived iron 
abundances. The star with the lowest k is HD180890. Lin- 
ear least square fits yield: 

a. k = -0.31 + 5.5 • 10~ 5 ■ T cS 

b. k = 0.23 - 0.054 -log g 

c. k = -0.014 + 0.073- [Fe/H]. 
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Fig. 6. Iron abundances, derived in this study, as a func- 
tion of effective temperatures. 



tential for the lines. 32 Fe I lines measured in most of the 
stars, marked in Table |[ were selected for this purpose. 
For all stars (excluding four stars with too few lines ob- 
served) the abundance of each iron line was plotted as a 
function of the excitation energy of the lower level, Xi > 
and a least square linear fit was made to the points of 
the form [Fe/H]= a + k ■ x„ Fig. § The mean value of 
the slope coefficients, k, is 0.00 (excluding the K dwarf 
stars), suggesting that the systematic errors in effective 
temperatures are small. 
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Next, these linear coefficients were plotted as func- 
tions of effective temperature, surface gravity and spectro- 
scopically derived iron abundance, Fig. ||. We have com- 
pared our results with those obtained by Edvardsson et 
al. (1993a) (their Fig. 9a-f) and find that the two studies 
span the same range of excitation energy slopes, k. There 
are small but systematic deviations from the excitation 
equilibrium, varying with effective temperature. A change 
in T cff of +100 K results in a Afc of +0.006 dex eV" 1 . The 
change of excitation energy slope with surface gravity esti- 
mates seems to reflect the variation of surface gravity with 
effective temperature (surface gravity increases as effective 
temperature decreases) . 

In many studies surface gravities are determined by 
requiring ionization equilibrium. This is typically made 
by changing the surface gravity until the iron abundances 
derived from Fe I and Fe II lines yield the same abundance. 
We have not, as discussed earlier, used this method. As 
we will see this has led us to discover what appears to be 
a case of significant overionization in K dwarf stars and 
an opposite smaller effect for the hotter stars (see Fig. Q). 

What would the effects be if we assumed ionization 
equilibrium, and adjusted the surface gravities accord- 
ingly? We can estimate changes in the stellar abundances 
from the results of Table || and Table [ll]. From this we find 
that half of the stars should have their surface gravities 
increased by 0.25 — 0.35 dex to achieve ionization equi- 
librium for iron. This means that the the iron abundance 
will change with ss —0.03. Abundances of other elements 
will change with similar amounts but with differing signs, 
which means that for some elements [X/Fe] will change 
by up to 0.1 dex and for others not at all. However, we 
note that the oxygen abundances are very sensitive to the 
surface gravity and may change by up to 0.2 dex. As a 
comparison we estimate the maximum error in the derived 
oxygen abundance caused by incorrectly set continua to 
be less than 0.1 dex. 

An adjustment downwards of the gravities by about 
0.3 dex would increase the conflict with the gravity esti- 
mates from the Cal 6162 A line wings. We consider such 
a revision less probable. 

5. 8. K dwarf stars - Overionization 

Our results admit a comparison for five elements (scan- 
dium, vanadium, chromium, iron and yttrium) of abun- 
dances derived from ions to abundances derived from the 
corresponding atom, as function of effective temperature, 
within a rather wide range of effective temperature. 

We find an apparent overionization as compared to 
expectations from LTE calculations for the five K dwarf 
stars in our study, Fig. |?]. Our results are at first sight 
unexpected, especially for the two K dwarf stars with iron 
abundances of ~ 0.3 dex relative to the solar iron abun- 
dance, as derived from Fei lines. In stars more metal-rich 
than the Sun the photoionizing non-local UV-flux is more 
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Fig. 7. Abundance ratios of chromium and iron as func- 
tions of [Fe/H]. [Crii/Cn], denotes [Cr/H] L TE as de- 
termined from Cm lines minus [Ct/ITjlte as deter- 
mined from Cri lines, and similarly for iron. Results for 
chromium are shown in panel a. and for iron in panel 
b.. x symbols denote the K dwarf stars and * symbols 
denote the stars from Barbuy & Grenon (1990). We ex- 
emplify, with HD61606A ([Fe/H] = -0.08)and HD103932 
([Fe/H] = 0.16), how the resulting loci of the K dwarf stars 
will be shifted in the diagram if the parameters of the stel- 
lar model atmospheres are changed; changes of T c g + 200K 
(solid lines), logg±0.4 dex (dashed lines) and [Fe/H] ±0.2 
dex (dotted lines) are shown. Filled triangles denote in- 
creased values of the respective parameters while filled 
squares denote decreased values. 



efficiently blocked than i n mo re metal-poor stars. How- 
ever, as discussed in Sect. 5.9, the stronger blocking may 
be more than compensated for by the increased tempera- 
ture gradient which may enhance the non-locality of the 
radiation field. 

We have carefully inspected the Cr II and Fe n lines in 
the K dwarf spectra and excluded all lines which may be 
subject to severe blends in these cool stars. For Fell we 
used the lines at 6456.39 and 6516.39 A and for Cm the 



12 



Sofia Feltzing & Bengt Gustafsson: Abundances in metal-rich stars 



Table 4. The table shows results of different derivations of the effective temperatures for the five K dwarf stars; in 
the second column as derived from the calibration by Olsen (1984), and which are used in this study; in the third 
spectral classes; in the fourth and fifth differences between effective temperature derived in this work and from MK 
class calibrations by Bell & Gustafsson (1989) and Johnson (1966); in the sixth and seventh B — V and effective 
temperatures derived from B — V using the calibration of Johnson (1966). The three last columns contain effective 
temperatures derived by other authors, Neff et al. (1995), Morell (1994) and Arribas & Martinez Roger (1989). 
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ID 


T e ff 




ATrff 


AT off 
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T e ff 


T e f£ 


HD32147 


4625 


K3V 


+200 


+100 


1.06 


4619 


4570 




4670+150 


HD61606A 


4833 


K2V 


+300 


+100 


0.96 


4863 








HD103932 


4510 


K5V 


-300 


-100 


1.12 


4473 








HD131977A 


4585 


K4V 


-100 





1.10 


4522 


4575 


4570 




HD136834 


4765 


K3V 


+100 





1.00 


4765 









lines at 5305.86 and 5310.69 A. In spite of using blend free 
lines the apparent overionization remains. The internal 
consistency between abundances derived from the two Fe II 
lines is very good and this is also the case for Cr II. 

For scandium, vanadium and yttrium blends remain a 
possible source of error, but the similarity with the trends 
for iron and chromium suggests a common cause of the 
apparent overionization for all these elements. 

The iron abundances derived from the atom show no 
obvious trend with effective temperature, see Fig. ^| 

A pro babl e reason for these effects is overionization 
(see Sect. 5J)). Before discussing this, however, we shall 
explore the possibility that errors in the temperature scale 
could also contribute significantly. 

5.8.1. Errors in effective temperatures for K dwarf stars 

The calibration of the photometry in Olsen (1984) is, for 
b — y > m .510, based on a sample of 15 K and M dwarf 
stars using stellar parameters from Cayrel de Strobel & 
Bentolila (1983). The K and M dwarf stars in Olsen's sam- 
ple span a small range in 8m\ and 8c\. This is reflected 
in the change in the calibration at b — y = m .510. For 
b—y> m .510 the calibration is degenerate in metallicity. 
Olsen (1984) quotes an error of ~100K for the effective 
temperature as derived from b — y. 

Additional photometry is scarce for our K dwarf stars. 
However, UBV photometry exists and we have checked 
our effective temperatures using the calibration of B — V 
by Johnson (1966), Table |]. These effective temperatures 
agree well with those obtained from uvby — (3 photometry. 
We note, however, that the increased blocking for metal- 
rich stars, as compared to the calibration stars used by 
Johnson (1966), may cause the present effective tempera- 
tures to be underestimated. We have also derived (crude) 
effective temperatures from calibrations of the MK classifi- 
cation. Comparing the effective temperatures derived from 
photometry with calibrations of effective temperatures as 



functions of spectral classification (Bell & Gustafsson 1989 
and Johnson 1966) we estimate T G ff-values that deviate as 
compared with our standard values as indicated in Table 
[|, columns 4 and 5. 

We have searched the literature for independent 
derivations of the effective temperature for these stars. 
Those found agree well with the photometric results, Table 
||. As described earlier we have also used the excitation en- 
ergy balance to check our effective temperatures. A change 
in effective temperature of +200K for HD32147 changed 
the slope of the least-square fit to the data points in the di- 
agram abundance-versus-excitation energy from +0.04 to 
+0.00. A change of -200K brought about a similar change 
but in the opposite direction. Thus, the excitation equilib- 
rium indicates that +200K is an acceptable change of the 
effective temperature for this star. The line-to-line scat- 
ter in derived abundances from Fei lines for this star is 
among the largest, ±0.13. The same changes in effective 
temperatures give similar values for HD61606A. Since the 
other three K dwarf stars are in the same effective tem- 
perature range and have similar surface gravities as the 
two stars discussed here, and since HD32147 is the most 
metal-rich and HD61606A is the most metal-poor of the 
five K dwarf stars, changes in the fundamental parameters 
of the remaining three stars will produce similar changes 
in abundances. 

To conclude, a change in effective temperature of +200 
to +400K may be allowed as judged from the excitation 
equilibrium. As is obvious from Fig. [?], a change of this 
order of magnitude would restore the LTE ionization bal- 
ance for iron and chromium. The large linc-to-linc scatter 
in derived abundances for lines with high excitation ener- 
gies makes attempts to derive effective temperatures from 
excitation equilibria very dependent on one or two points 
in the lower end of the excitation energy range, spanned 
by the lines as illustrated in Fig. ^. Effective temperatures 
derived in other studies and from B — V colours deviate 
by less than this from our values. 
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As already noted the surface gravities for the K dwarf 
stars seem rather well determined, see Table § We also 
note that our \ogg value for HD131977A agrees well with 
that given by Morell (1994). We conclude that realistic 
errors in log g are not enough to account for the departure 
from ionization equilibrium. 

To conclude, we cannot from our analysis exclude that 
the apparent pattern of overionization, at least partially, 
is due to a temperature scale that is several hundred K too 
low. However, our analysis, together with evidence from 
other studies, suggest that deviations from LTE is a more 
plausible cause for the effects. 

5.9. Non-LTE 

No detailed study has been devoted to the non-LTE effects 
on abundance determinations for metal-rich dwarf stars, 
cooler than the Sun. A general result of the available stud- 
ies for solar-type stars is, however, that several different 
effects are at play and may counteract each other, and this 
makes all extrapolation to the present study of metal-rich 
dwarf stars from studies of other types of stars or studies 
of other elements questionable. 

Among the significant effects are (cf. Bruls et al. 1992) 
resonance line scattering, photon suction, ultraviolet ove- 
rionization, (infra)red over-recombination and optical line 
pumping. The resonance-line scattering effects, in which 
photon losses cause the source functions of resonance lines 
to drop far below the Planck function at depths greater 
than those where the line optical depth is unity, may lead 
to severe overestimates of abundances - e.g., Carlsson et 
al. (1994) find that in Li-rich cool stars the Li abundance 
may well be overestimated by a factor of 3 as a result of 
this. 

Photon suction may, for metal-rich cases in particu- 
lar, lead to overpopulation of, e.g., the ground state and 
thus inhibit the effects of overionization. This is the result 
of a compensation of photon losses in the upper photo- 
sphere in resonance lines, as well as in connected ladders 
of transitions, by a downward population replenishment 
flow from the continuum reservoir. It is of great signif- 
icance for atoms with a majority of corresponding ions 
and with pronounced cascade ladders. For complex atoms 
is should be of greatest significance for the high-lying levels 
that thus can compensate population depletion processes 
at lower excitation energy, e.g. caused by overionization. 

The ultraviolet overionization has been a major worry 
in analyses of late type stellar spectra for two decades. 
It arises because the mean intensity J„ drops below the 
Planck function J3„ in the line-forming regions of the at- 
mosphere on the blue side of the spectrum peak. Overi- 
onization is known to occur for Fe I in the Sun from lev- 
els a few electron volts below the continuum (see Rutten 
1988 and references therein) and may more or less effect 
other metals as well (see, e.g., Baumueller & Gehren 1996, 
Bruls 1993). Overionization was suggested by Auman & 



Woodrow (1975) to be significant for a number of ele- 
ments with lower ionization energies in cool stars. Major 
problems in modeling it are, however, the difficulties in 
predicting the ultraviolet flux of late-type stars with the 
crowding of spectral lines and the possible existence of an 
"unknown opacity" (cf., e.g. Gustafsson 1995), as well as 
the dependence of the results of the notoriously uncertain 
collision cross sections, e.g. for collisions with H atoms 
(cf. Steenbock & Holweger 1984). Empirically, Mackle et 
al. (1975) and Ruland et al. (1980), found a tendency for 
the low-excitation lines (x ~ 2 eV) of Fe I and Ti I in K gi- 
ant star spectra to give systematically lower abundances 
than the high-excitation lines (x ~ 4.5 eV). The abun- 
dance difference was typically found to be 0.15 dex. Sub- 
sequently, Steenbock (1985) succeeded in reproducing this 
result with statistical-equilibrium calculations. He found 
the effect to mainly reflect an overionization in upper lay- 
ers (notably T5000 < — 2) of the atmospheres, where the 
low excitation lines are formed. The effect is much smaller 
in the solar spectrum, leading to systematic errors in a 
differential analysis where red giants are compared with 
the Sun. 

For metal-rich stars, the blocking by the crowd of spec- 
tral lines in the ultraviolet could be expected to - at first 
sight - strongly reduce the overionization effect, but this 
may be compensated for by a steeper temperature gradi- 
ent in their atmospheres as a result of line blanketing ef- 
fects, as in the case of Lii, Carlsson et al. (1994), or of Cai, 
Drake (1991). The latter study is particularly instructive 
for judging the results of the present investigation. Drake 
finds that for G and K-type stars the overionization effects 
on Ca 1 abundances increase with decreasing effective tem- 
perature, with increasing acceleration of gravity and with 
increasing metallicity. At least the last two results may 
seem contrary to intuition. They reflect the significance 
of H 1 absorption shortwards of the Balmer discontinuity, 
which blocks more of the ionizing UV flux for the giants 
than for the dwarfs, and the afore-mentioned effects of 
metal-line blanketing on the temperature structure. For 
the K dwarf stars the effects on Ca 1 abundances may, ac- 
cording to Drake's results, well result in an underestimate 
by a factor of two or more if LTE is assumed. 

Over-recombination is important for photoionization 
transitions from levels close to the continuum (i.e.. in the 
infrared), since for them the angle-averaged J„ may drop 
below the local Planck function deep in the photosphere. 
This may produce net recombinations, and overpopulation 
of the upper levels. 

Optical (ultraviolet) pumping occurs in strong lines, 
e.g. the resonance lines and is analogous to overioniza- 
tion in that it is driven by J„ > B u . It is important, not 
the least when it occurs in ultraviolet resonance lines and 
excites the atoms to states which may be much more eas- 
ily photoionization due to a much richer radiation fields 
available at longer wavelengths, as was early suggested by 
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Aumann & Woodrow (1975). This is most important for 
trace elements, and for metal-poor stars. 

The complex interplay between these different mecha- 
nisms affects most levels of the atom, at great atmospheric 
depths, for atoms where the strong lines get efficiently op- 
tically thin in the photosphere, i.e. for relatively rare ele- 
ments like the alkalis, while for the more abundant atoms 
like Fe and Mg the stronger, e.g. resonance, transitions are 
in detailed balance through most of the photosphere. For 
these, the relatively simple overionization phenomenon is 
probably dominating, except for transitions very close to 
the continuum for which e.g. photon suction may be sig- 
nificant. 

In a recent study Gratton et al. (in prep.) have used 
detailed statistical-equilibrium calculations to explore the 
departures from LTE for solar-type dwarfs as well as for 
red giants of different metallicities, and their effects on 
abundance determinations for O, Na, Mg, and Fe. These 
authors find relatively small effects for stars cooler than 
the Sun for O I - LTE abundances from the IR triplet 
lines should be corrected downwards by less than 0.1 dex 
for stars with T e g ^ 6000 K. For Nal the subordinate 
lines are weakened by overionization and cascade by about 
0.1 dex for the solar- type dwarfs. The dominating effect 
for Mg I is overionization, and the non-LTE abundance 
corrections are thus generally positive. Typically the cor- 
rections are 0.1 dex in the dwarf stars. For Fei, where 
again overionization is dominating the abundance correc- 
tions ~ 0.1 dex. Most of these effects are found to be 
greater for T c s > 6000 K. It should be noted, however, 
that models for metal-rich dwarfs with T c g < 5000 K were 
not included in this study. 

Summing up the discussion of non-LTE we conclude 
that the effects on abundances are expected to be mainly 
due to overionization for most of the elements. For the 
alkali atoms, as well as for the rare earths, more complex 
effects may also be significant. Typically, errors of about 
0.1 dex may be expected in the differential results but the 
complexity of the interplay between different effects, and 
in particular the results obtained by Drake (1991) for Cai, 
suggest that greater effects may be present, in particular 
for the metal-rich K dwarf stars. 

5.10. Collecting errors 

We have shown that errors in fundamental parameters give 
errors in mean resulting abundances of less than 0.1 dex. 
For elemental abundances derived from several lines this 
may be the dominating error, while for abundances de- 
rived from one single line errors due to blends and fitting 
of continua may be the main contributors to the overall 
error. Deviations from LTE in the excitation and ioniza- 
tion balance may also be of importance, probably more so 
for abundances based on few lines, in particular for the K 
dwarf stars. We collect our best estimates of errors due to 
different sources in Table & 



Table 5. The effects of error sources explored in this work 
on estimates of abundances relative to the Sun. 



Source of error Error in resulting relative abundance 

Measurement of W\ negligible 

Continuum fitting < 0.09 dex, usually 0.05 dex 

T cS ,\ogg < 0.1 dex 

Non-LTE effects 0.1 - 0.2 dex? 

Oscillator strengths <0.1 dex 



Table 6. Comparison of iron abundances between our 
work and the abundances quoted in the catalogue by 
Cayrel de Strobel et al. (1997). The second column gives 
our results and the third the mean, and the spread, of the 
iron abundances given by Cayrel de Strobel et al. In the 
fifth give the number of derivations used. A straight mean 
has been taken to represent the mean abundance from the 
catalogue. 



ID 


[Fe/H] 


<[Fe/H]> 


±s 


# 


HD30562 


0.19±0.09 


0.14 


0.0 


2 


HD32147 


0.22±0.13 


0.02 


0.0 


2 


HD67228 


0.16±0.08 


0.05 




1 


HD131977A 


0.00±0.12 


0.01 




1 


HD 144585 


0.27±0.05 


0.23 




1 


HD182572 


0.42±0.05 


0.32 


0.14 


7 


HD186427 


0.12±0.05 


0.06 


0.04 


5 



5.11. Comparison of results for stars in common with 
other studies 

The majority of our stars have not been studied before 
through spectroscopic abundance analysis. 

For those of our stars (HD30562, HD32147, HD67228, 
HD1319777, HD182572, HD186427) that are in the cat- 
alogue by Cayrel de Strobel et al. (1997) the agreement 
between iron abundances derived in this study and those 
listed in the catalogue is good, cf. Table |[ 

HD32147 has been given much attention in the discus- 
sion of Super Metal Rich (SMR) stars. SMR stars have 
been defined as stars with [Fe/H] < 0.2 dex (for a discus- 
sion and references on SMR stars see Taylor 1996). Low 
resolution work and photometric determinations of [Fe/H] 
have been carried out for this stars, but this is, to our 
knowledge, the first high dispersion analyses of the star. 
Our [Fe/H] of 0.28 dex implies that this star is really an 
SMR star. 

Our results for HD182572 are compared with the re- 
sults of the detailed analysis by McWilliam (1990) in Table 
. We note that considerable discrepancies remain even af- 
ter correcting for the difference in effective temperatures. 

We have analysed three stars previously studied by 
Edvardsson et al. (1993a). The results are compared in 
Table |. For HD30562 and HD144585 the results agree to 
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Table 7. Comparison of abundances, [X/H], de- 
rived for HD182572 (HR7373) in our study and by 
McWilliam (1990). McWilliam uses (T cff /log 5 /[Fe/H]/&) 
= (5380/3.92/0.15/1.9) and we (5739/3.83/0.42/1.9). The 
results by McWilliam have been scaled to the same solar 
abundances as we use, Table |2|, (this is most important 
for iron). In the last column our values are scaled to the 
effective temperature used by McWilliam. 



Mc William This work 

5739 K 5380 K 



Fei 


0.31 


0.42 


0.24 


Si i 


0.28 


0.51 


0.47 


Cai 


-0.11 


0.42 


0.21 


Sen 


0.14 


0.36 


0.36 


Til 


-0.02 


0.50 


0.21 


Vi 


-0.02 


0.44 


0.12 


Co i 


0.18 


0.58 


0.29 


Nil 


0.00 


0.46 


0.24 


Eull 


0.18 


0.13 


0.13 



Table 8. The first column for of each star contains 
the result of this work, the second the difference be- 
tween the two studies, Diff = this work - Edvardsson 
et al. (1993a). Stellar parameters used by Edvardsson et 
al. are for HD30563 (T cff /log ff /[Fe/H])=(5886/3.98/0.17), 
for HD67228 (5779/4.20/0.04) and for HD144585 
(5831/4.03/0.23). The values of the parameters used in 
this study may be found in Table [|. 



HD30562 HD67228 HD144585 

Diff. Diff. Diff. 



[Oi] 






0.21 








Nai 


0.21 


+0.01 


0.23 


+0.06 


0.36 


±0.00 


Mgi 


0.33 


-0.01 


0.23 


+0.08 






All 


0.25 


-0.07 


0.22 


+0.03 






Si i 


0.21 


-0.02 


0.27 


+0.11 


0.24 


-0.03 


Cai 


0.17 


-0.01 


0.15 


+0.12 


0.24 


-0.03 


Til 


0.13 


-0.02 


0.09 


+0.02 


0.29 


±0.02 


Fei 


0.19 


+0.05 


0.16 


+0.12 


0.27 


±0.04 


Fen 


0.16 


+0.10 


0.22 


+0.06 


0.11 


-0.05 


Nil 


0.17 


±0.00 


0.14 


-0.01 


0.32 


±0.06 


Yn 


0.09 




0.12 


±0.00 


0.06 


±0.06 


Ndn 






-0.10 




-0.14 


-0.06 



Table 9. Comparison of derived abundances, [X/H], be- 
tween Friel et al. (1993) and this work for HD186427, 
16 Cyg B. Friel et al. use the following stellar pa- 
rameters (T off /log 5 /[Fe/H])=(5770/4.30/0.05) and we use 
(5773/4.17/0.12). 





Nai 


All 


Si i 


Cai 


Friel et al. 
This work 


0.07 

0.13±0.04 


0.12 

0.12±0.04 


0.06±0.03 
0.10±0.04 


0.07±0.05 
0.02±0.09 




Til 


Fei 


Fen 


Nil 




0.10±0.04 
0.07±0.10 


0.05±0.03 
0.05±0.05 


0.02±0.04 
-0.06±0.04 


0.05±0.03 
0.04±0.10 



Table 10. Comparison between results of our anal- 
yses for HD30592 and HD67228, with the results of 
Tomkin et al. (1997). The second and fourth column 
give the results of this work, while the third and fifth 
give the difference between this work and that of Tomkin 
et al. Stellar fundamental parameters of Tomkin et al. 
were (T cff / log ff /[Me/H]/£ t ) = (5930/4.20/0.21/1.55) and 
(5835/4.10/0.09/1.75) for HD30592 and HD67228, respec- 
tively. Our parameters are given in Table @. 





HD30562 


HD67228 






Diff. 




Diff. 


[Oi] 




0.21 


0.01 




Nai 


0.21 


-0.02 


0.23 


±0.03 


Mgi 


0.33 


0.12 


0.23 


0.09 


All 


0.25 


-0.01 


0.22 


0.03 


Si i 


0.21 


-0.04 


0.27 


0.12 


Cai 


0.17 


-0.06 


0.15 


0.08 


Sen 


0.31 


0.06 


0.33 


0.16 


Til 


0.13 


-0.12 


0.09 


-0.02 


Vn 


0.10 


-0.16 


0.12 


0.01 


Cri 


0.18 


-0.02 






Cm 


0.17 


-0.05 


0.18 


0.02 


Fei 


0.19 


-0.02 


0.16 


0.07 


Fen 


0.16 


-0.04 


0.22 


0.10 


Nil 


0.17 


-0.08 


0.14 


0.01 


YIi 


0.09 


-0.12 


0.12 


0.06 


Zn 


0.65 


0.39 






Eull 


0.14 


-0.17 







within the errors quoted. The results for HD67228 show 
larger discrepancies than the other stars for magnesium, 
silicon, calcium and iron. A higher microturbulence pa- 
rameter, as used in Edvardsson et al. (1993a), would de- 
crease our results by 0.01 or 0.02 dex (ionized iron by 0.04 
dex). Edvardsson et al. (1993a) find that a change in T e g 
by ±100 K gives an iron abundance 0.06 higher for their 
stars. Thus, the increase of 52 K needed to transform the 
results of Edvardsson et al. (1993a) to our temperature 



scale means an increase of their iron abundance by 0.03 
dex. Also the silicon and calcium abundances are affected 
in the same way as iron while all the other abundances 
remain as before, within the errors, in the two studies. 
The log<?/-values agree well (0.06 difference, we have the 
higher value) . The spectrum we have obtained for this star 
is of high quality (S/N - 200). We do not find the dis- 
crepancy between the two studies alarming. HD67228 has 
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also been studied by Andersen et al. (1984) in a study on 
lithium isotope ratios in F and G dwarf stars. They derive 
an [Fe/H] of 0.05 from spectral lines using a model with 
(T cff /log 5 /[Fe/H])=(5850/4.2/0.05). 

Two stars from our sample, HD30562 and HD67228, 
were recently analysed in detail by Tomkin et al. (1997) on 
the basis of different spectra; however, obtained with the 
same instrument and analysed independently with model 
atmospheres computed with the same computer program. 
The results are compared with those of our analyses in 
Table [l^. In view of the errors in these analyses we find 
the agreement satisfactory. 

HD186427 (16 Cyg B) have been extensively studied, 
in particular in connection with searches for solar twins. 
A recent spectroscopic study has been performed by Friel 
et al. (1993). The results are in good agreement, Table ^|. 
However, we find lower iron abundances derived from Fen 
lines than Friel et al. (1993) do. This difference is probably 
mainly due to the different surface gravities used. 

To conclude, we find that, for those few stars in our 
programme in common with other studies, abundances are 
rather well reproduced. This gives confidence when we now 
apply our results to the exploration of the chemical evo- 
lution of the Galaxy. 

6. Abundance results 

Our abundance results are presented in Table |ll| and rel- 
ative abundances are plotted in Figs. [§[ |§ andg. Below, 
we shall compare our results with model calculations for 
the galactic chemical evolution. Sometimes, we shall also 
quote results of other studies, in particular that of Ed- 
vardsson et al. (1993a) (which is grossly compatible in 
terms of calibration and methods with the present one), 
but also with others, in particular for Population II dwarf 
stars. 

Comparison of stellar abundance results from differ- 
ent studies may not always be straightforward. There are 
many inconsistencies that may confuse the interpretation 
of such comparisons. The usage of different lines for the 
abundance analysis, different stellar parameters, different 
model atmospheres, and, if the study is differential to the 
Sun, different solar model atmospheres; all of these may 
lead to offsets and may cause the compiled data to show 
trends which are unreal. Nevertheless, below we shall com- 
pare our abundances with results from other studies to put 
our results into a broader picture of the galactic chemical 
evolution. 

6.1. Iron 

In observational studies of galactic chemical evolution iron 
is often used as the reference element. The reason for this 
is twofold; iron is believed, but this is debated, to be a fair 
chronometer for the nucleosynthesis in the Galaxy, and 
the spectra of dwarf stars show many iron lines, easy to 



measure. The evolutionary picture for iron is complicated 
by the fact that iron is produced in both core collapse 
and type la supernovae. From this point of view oxygen, 
which is only produced in core-collapse supernovae, may 
be preferable as reference element. However, as we will 



discuss in Sect. 3.2, oxygen abundances are not trivial to 



derive. We will therefore conform with common practise 
and use iron as reference element. 

Our resulting iron abundances are well determined 
with a line-to-line scatter in [Fe/H] of typically 0.09 dex 
and a formal error in the mean iron abundance of typically 
less than 0.02 dex for each star. We find that [Fe/H] does 
not vary with Vlsr, i.e. R m , Fig. |[ This is an important 
observation when we consider other elemental abundances 



relative to Fe later, cf. Sect. S.3, and a first indication that 
the mixing of gas over the Ai? m spanned has been quite 
efficient. The five K dwarf stars show a similar behaviour 
as the rest of the sample. 

6.2. Oxygen 

Oxygen is the third most abundant element in stars and 
therefore plays a significant role for stellar opacities and 
energy generation. Therefore, determination of stellar ages 
depends strongly on the assumed initial oxygen abun- 
dance in the star, see e.g. VandenBerg (1992). Oxygen 
abundances also affect the determination of time-scales in 
the galactic chemical evolution and star-formation rates. 
Thus, it is important to know the amount of oxygen 
throughout the history of the Galaxy. 

We have studied three oxygen criteria; the forbidden 
line at 6300A, the 6158A line and the 7774A triplet. These 
criteria are commonly used; however, discrepancies be- 
tween the abundances derived from the different crite- 
ria for the Sun, as well as for other late-type stars, have 
prevailed in spite of much work (cf. e.g. Eriksson & Toft 
(1979) and Kiselman (1993) and references therein). Non- 
LTE and granulation are two proposed sources of the dis- 
crepancy between abundances derived from the [O i] line 
and from the triplet lines. 

The formation of the [O i] line is expected not to be 
subject to departures from LTE. The lower level of the 
transition is the ground level of the atom and the major- 
ity of the oxygen should be found in the ground state of 
the atom under solar photospheric conditions. There are, 
however, suspicions that the analysis of the line might be 
subject to systematic errors due to the adoption of plane 
parallel stellar atmosphere models since granulation ef- 
fects are not taken into account in these models, Kiselman 
& Nordhmd (1995). 

Edvardsson et al. (1993a) found a correlation between 
oxygen abundances derived from the [O i] line and the 
abundances derived from the triplet lines as well as a cor- 
relation with the abundance derived from the 6158A line. 
The se relations are shown, together with our data, in Fig. 
.21 We do not find a clear correlation between abundances 
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Table 11. Elemental abundances our programme stars. For each star and each ion the derived elemental abundance, 
the line-to-line scatter (the error in the mean = line-to-line scatter /valines) an( ^ the number of lines used in the 
analysis are given. The table is continued on the following pages. This table is only published electronically. 
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Fig. 8. Our abundances relative to Fe. x symbols denote the five K dwarf stars while * symbols denote the stars from 
Barbuy & Grenon (1990). 



derived from [O i] and the abundances derived from the 
7774A triplet and 6158A lines for our stars. This circum- 
stance, and the fact that the [0 1] line is not expected to 
be affected by departures from LTE, lead us to use only 
the abundance derived from the [Oi] line in our analysis. 

When comparing our results with those of Nissen & 
Edvardsson (1992) we find in general a higher oxygen 
abundance. This discrepancy can be understood: Nissen 
& Edvardsson (1992) have used an oscillator strength of 
-9.75 (Lambert 1978) while ours is -9.84. This means that 
our abundances should be scaled down by 0.09 dex to be 
put on the same scale as the Nissen & Edvardsson (1992) 
abundances. 

Three of our stars, HD37986, HD77338 and HD87007, 
have previously been studied by Barbuy & Grenon (1990). 



These authors derive oxygen abundances and metallici- 
ties for a group of 11 dwarf stars. The stars were selected 
on the basis of their kinematics and claimed to represent 
the "local bulge population" . The stars fell clearly above 
the oxygen trend expected from simple models of galactic 
chemical evolution (their Fig. 1). The results were inter- 
preted as possible evidence for a rapid, and probably early, 
enrichment of the gas in the galactic Bulge. For most of 
the 11 stars in their study, Barbuy & Grenon (1990) de- 
rived the same abundance from the forbidden line as from 
the triplet lines. This is not the cas e for the majority of 
the stars in our study, see Fig. 6.2. As shown in Table 
for the three stars in common with Barbuy & Grenon 
(1990), we find that our iron abundances are lower by 0.2 
dex as compared to their results, while the oxygen abun- 
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Fig. 9. Our abundances relative to Vlsr- x symbols denote the five K dwarf stars. 
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Fig. 10. Iron abundances relative 
hydrogen vs Vlsr- x signs denote 
the five K dwarf stars and HD87007 
is denoted by a * sign. 



dances derived from the triplet lines stay the same relative 
to iron. The largest discrepancy in derived oxygen abun- 
dance is found for the forbidden line in HD87007. Unfor- 
tunately, the error in our determination of this abundance 
is rather large. The spectrum is one of our poorer, with 
a S/N of only ~ 80. The estimated error in [O/H] from 
noise might then be as large as 0.5 dex. Thus, it is possible 
that the high abundance derived by us from the forbid- 
den line for HD87007 is due to errors. For the two other 



stars the forbidden oxygen line was, unfortunately, heavily 
obscured by telluric lines and could not be used. 

We only have access to velocity data for one of the 
stars, HD87007 (Barbuy, private communication). It has 
a ^lsr velocity of —42.5 km/s which means that it would 
not satisfy the velocity criteria of our high-velocity sam- 
ple designed to represent the inner part of the disk. The 
high ZTlsr velocity (radial velocity in the galactic plane of 
symmetry) of the star, however, gives a total spatial ve- 
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Fig. 11. Oxygen abundances de- 
rived from different abundance cri- 
teria are compared, panel a., b. and 
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A oxygen abundance represents the 
mean of the results for the three 
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for different velocities. The relation 
between Vlsr an d R m is taken from 
Fig. 1 in Edvardsson et al. (1993b). 
The star with the lowest oxygen to 
iron abundance ratio is HD110010. 



locity which fulfils the requirement of membership in our 
high velocity sample. Castro et al. (1997) have studied 9 
other high velocity stars from the work by Grenon (1989). 
Their results support our results. It is, however, difficult 
to make a clear comparison since velocity data for these 
stars have not been published. 



In Fig. 6.2 we plot [O/Fe] vs Vlsr for our stars as well 
as the, generally, more metal-poor disk stars in the sample 
of Nissen & Edvardsson (1992). For the combined sample 
one may trace a tendency for [O/Fe] to decrease with in- 
creasing VlsRi but this is not shown by our sample alone. 
From Fig. [l^ we conclude that the oxygen abundance 
in general keeps declining relative to the iron abundance 
also for [Fe/H] > 0.1 dex. This is not inconsistent with 
what Nissen & Edvardsson (1992) found in this metallic- 
ity range. 

Oxygen is produced in massive stars exploding as su- 
pernovae of types II, lb and Ic, Woosley & Weaver (1995) 
and Thielemann et al. (1996). Therefore, oxygen is ex- 
pected to rapidly build up at early times in the Galaxy or 
in any region which has experienced substantial star for- 



Table 12. Stellar parameters and oxygen abundances for 
the stars in common between this work and that by Bar- 
buy & Grenon (1990). For each star we give our iron 
abundandance as derived from spectral lines, [Fe/H], the 
metallicity as quoted by Barbuy & Grenon, [M/H] and 
the oxygen abundances derived from the different crite- 
ria as indicated. For each criterion we give our resluts on 
the first line and the results by Barbuy & Grenon on the 
second. [0/Fe]6iss was not measured by them. 



ID 


HD37986 


HD77338 


HD87007 


[Fe/H] 


0.27 


0.22 


0.27 


[M/H] 


0.47 


0.45 


0.43 


[0/Fe] 63 oo 






0.54 




0.15 


0.20 


0.00 


[O/Fc] 7 770 


0.15 


0.29 


0.27 




0.23 


0.20 


0.20 


[0/Fe] 6158 


0.39 


-0.04 


0.27 
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Fig. 13. Oxygen abundances from this work, o symbols, 
and from Nissen & Edvardsson (1992), • symbols, x sym- 
bols denote our five K dwarf stars, * symbol HD87007 and 
the filled boxes data from Castro et al. (1997). The star 
with lowest [O/Fe] is HD1 10010. 

mation "lately" . [O /Fe] starts to decline once the iron pro- 
ducing supernovae start contributing more significantly to 
the enrichment of the interstellar gas. This decline starts 
at [Fe/H]~ —1 dex in the galactic disk, or at even lower 
metallicities. 

In Fig. [l3| we compare our data with three different 
theoretical models of the galactic chemical evolution. The 
model by Matteucci & Frangois (1989) clearly shows the 
envisaged decline in [O/Fe] after [Fe/H] = -1.0 dex, i.e. 
when supernovae type la start to contribute to the enrich- 
ment of the gas. This can be compared with a more recent 
model, taking the effects of metallicity dependent super- 
nova yields into account, by Prantzos & Aubert (1995). 
The difference between the two models by Prantzos & 
Aubert (1995) is not large for lower metallicities but from 
around solar metallicity there is an increasing discrepancy 
between their two models. If this suggested trend contin- 
ues to higher [Fe/H] the model using metallicity depen- 
dent yields would be favoured by our data. The model by 
Pagel & Tautvaislene (1995) is a simple analytic model, 
which assumes supernovae type la to give their yields af- 
ter a fixed time delay. In spite of its simplicity it fits the 
data remarkably well and may suggest that the basic un- 
derstanding of the processes involved is correct. 

Tsujimoto et al. (1995) studied the abundance gradi- 
ents in the galactic disk by means of a viscous disk model 
of galactic chemical evolution. As can be seen in their Fig. 
8 differences in oxygen abundances as a function of radial 
distance from the galactic centre are predicted to be so 
small that it would hardly be possible to resolve this in 
a study as ours where such a small part of the radius is 



suggested trend is continued when the models are further 
evolved, it would be difficult to reconcile them with our 
data, and with the results by Castro et al. (1997). 

6.3. Sodium 

One of the aims with our study was to determine whether 
or not the upturn of [Na/Fe] vs. [Fe/H] found by Edvards- 
son et al. (1993a) is real and if so, if it is an effect of a 
mixture of stars from different populations, cf Figs. [|, || [l4| 
and 6.3. We confirm that the upturn is real. However, the 



spanned. This is compatible with our results in Fig. 6.2. 
We note that their model predicts (their Figs. 7 and 8) 
[O/Fe] vs [Fe/H] to flatten out at solar metallicities. If this 



upturn is less steep in our study than what one could trace 
from the scattered diagram of Edvardsson et al. (1993a). 

In Fig. |6.3| we have divided our sample into stars rep- 
resenting the disk interior to the solar orbit and the solar 
orbit and made linear least-square fits to the data. In panel 
a. we show all 47 stars and a least-square fit with the K 
dwarf stars and the stars from Barbuy & Grenon (1990) 
excluded. In panel b. we show the stars that represent 
the disk interior to the solar orbit. There is no apprecia- 
ble difference found between this sample and the whole 
sample. In panel c. and d. we have defined the solar or- 
bit sample in two different ways. In c. it contains all stars 
with Vlsr < 45 km/s and in d. with Vlsr < 30 km/s. 
The sample in panel c. also has a behaviour indistinguible 
from that of panel a. However, the Vlsr < 30 km/s seem 
to show a somewhat steeper trend. 

An interesting question is now what difference in 
[Na/Fe] one would have reason to expect for stars 
formed at different R m . Using the data of Edvardsson 
et al. (1993a) we estimate that the minimum value of 
d[Na/Fe]/dR m is about 0.05 dex/kpc. Extrapolating this 
to the metal-rich stars and to i? m ~ 6 kpc one finds at 
the most a difference of 0.1 dex between our two samples. 
Implicit in this assumption is then that the population 
of metal-rich stars at 6 kpc, which is very sparsely repre- 
sented in the sample of Edvardsson et al. (1993a), is not 
qualitatively different from that in the solar neighbour- 
hood. The results obtained in the present study supports 
this and indicate that the difference in [Na/Fe] is, in fact, 
at the most 0.05 dex. 

In this connection, one should also note that orbital 
diffusion may well mask the possible differences between 
stars formed at 6 kpc and 8 kpc. E.g., stars of solar age in 
an orbit with R m — 8.5 kpc may have migrated from an 
orbit with R m = 6 - 7 kpc (cf. Wielen et al. 1996). The 
mixture of stars with different original orbits, in combi- 
nation with a radial galactic gradient of [Fe/H], was pro- 
posed by Wielen et al. (1996) to explain the unexpectedly 
high scatter in [Fe/H] of 0.2 dex for solar type stars, with 
similar age and similar present R m , found by Edvardsson 
et al. (1993a). From the work of Wielen et al. (1996) we 
estimate that two samples of stars with < R m >= 8.5 and 
6.5 kpc, respectively, would then be mixed by orbital dif- 
fusion so much that the population effects in abundances 
only would show up to about half the expected size as com- 
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the five K dwarf stars which are denoted by x symbols and the three stars from Barbuy & Grenon (1990) which are 
denoted by * symbols. Results from Edvardsson et al. are denoted by • symbols. 



pared with the situation if orbital diffusion is not present. 
Although the reason for the great inhomogeneities in the 
gravitational potential, needed to account for the orbital 
diffusion of this magnitude, is not known, we conclude 
that the effects looked for by dividing the total sample of 
stars according to the velocity criteria used here, might be 
diminished considerably by this phenomenon. 

6.4- Aluminium 

Aluminium is produced in heavy single stars, with M > 
8M Q . The scatter, as well as the internal line-to-line scat- 
ter, in our data is considerably smaller than in Edvardsson 
et al. (1993a). The upturn in [Al/Fe] vs [Fc/H] indicated in 
their data is not obviously present in ours, see Fig. [l4|. The 
two studies use different lines for the abundance analysis. 
The small scatter in [Al/Fe] for [Fe/H] > 0.0 dex is also 
evident in Morell (1994). The lines used by Morell (1994) 
and us, 6696.03A and 6698.66A, are situated in a part of 



the stellar spectrum which is clean. Thus, we expect no 
problems with continuum fitting and blends and errors 
arising from measurements of the line strengths should be 
negligible. The trend in Fig. [w] and its similarity with, 
e.g., that of Ca (Fig. |l8|) suggests a similar origin in core- 
collapse supernovae. 

6.5. Magnesium 

Due to the large line-to- line scatter in our magnesium data 
it is not possible to determine here if the large scatter 
found by Edvardsson et al. (1993a) is real or not. We note, 
however, that we get roughly the same amount of scatter 



as they do, Fig. 16. There is no evidence in our data for a 



correlation between kinematics and [Mg/Fe] ratios of the 
stars, Fig. ^. 

Magnesium is, like oxygen, usually assumed to be 
formed only in core-collapse supernovae through hydro- 
static carbon burning. Timmes et al. (1995) are not so 
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Fig. 15. In this figure we explore 
the different velocity samples as 
indicators of differences in galac- 
tic chemical evolution in adjacent 
parts of the galactic disk. The five 
K dwarf stars and the stars from 
Barbuy & Grenon (1990) are ex- 
cluded from this discussion. In panel 
a. we show [Na/Fe] for all stars. 
Error-bars are smaller than the sym- 
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Vlsr < -50 km/s and/or Q LS r > 
60 km/s shown. In panel c. stars 
with Qlsr < 45 km/s and in panel 
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linear least square fit to the data 
and the dashed lines show least dis- 
tance fits to the same data. Stars 
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Fig. 16. Magnesium vs iron abundances from several studies, as given in the figure, as well as model calculations of 
the galactic chemical evolution from Pagel & Tautvaisiene (1995) (dashed line) and Timmes et al. (1995) (solid line), 
x signs denote the K dwarf stars in our sample. Indicated error bars refer to the error in the mean. 



successful in describing the over-all evolution of magne- 
sium abundances in the Galaxy. This may suggest that 
the production of Mg is not fully understood at present. 
The simple-minded model with the delayed yield formal- 
ism is more successful in this respect, Pagel & Tautvaisiene 
(1995). 



6.6. Silicon 

We only use two lines to determine silicon abundances 
while Edvardsson et al. (1993a) used eight lines. We find 
a much larger scatter in our data than they do, Fig. [l4| 
It has not been possible, from our data, to determine the 
origin of neither the line-to-line scatter nor the star-to- 
star scatter. When inspecting the lines one by one no line 
seems to stand out in terms of derived abundances. Nor do 



Sofia Feltzing & Bengt Gustafsson: Abundances in metal-rich stars 



23 




0.6 



0.4 
0.2 
0.0 

-0.2 
-0.4 



1 

b. 




□ 

A 


X 


X 




- * 

X 


□ 


□ 


4500 


5000 


5500 6000 



eff 



Fig. 18. Calcium abundances for our stars, except the K 
dwarf stars, o symbols. The K dwarf stars are denoted by 
x symbols. Abundances for solar-type stars from Abia et 
al. (1988), □ symbols, and Gratton & Sneden (1987), A 
symbols, and the stars from Barbuy & Grenon (1990), * 
symbols. 



we find any signs of the scatter to be an effect of different 
stellar populations mixing, see e.g. Fig. |9[ 



6. 7. Calcium 

Comparison of our calcium abundances with those of Ed- 
vardsson et al. (1993a), in Fig. confirms their finding 
that the [Ca/Fe] flattens out towards higher metallicities. 
Also, there is no obvious difference between stars with dif- 
ferent galacto-centric mean distances, Fig. 

As is evident from Figs. || and [l^ the K dwarf stars ex- 
hibit a behaviour which is very different from the rest of 
the sample. They seem to have a mean calcium relative to 
iron abundance ~ 0.2 dex lower than the mean abundance 
for the rest of the stars. When plotting [Ca/Fe] as a func- 
tion of effective temperature we see that the stars with low 
Ca abundances have the lowest effective temperatures, cf. 
Fig. [1§ 

In a high resolution, high signal-to-noise abundance 
study of dwarfs and giants in the disk by Abia et al. (1988), 
there are six dwarf stars in the same metal and effective 
temperature range as our K dwarf stars. In Fig. 4a in 
Abia et al. (1988) the same features as in Fig. [1^ can be 
seen, i.e. cool, metal-rich stars show up as underabundant 
in calcium. In another study, Gratton & Sneden (1987), 
of light elements in field disk and halo stars we also find 
support for such a behaviour of calcium in cool dwarf stars. 



The key to these low [Ca/Fe] measures may lie in ove- 
rionization. Drake (1991) performed non-LTE calculations 
for calcium for a range of stellar parameters. He showed 
that the difference between an abundance derived under 
non-LTE and LTE conditions varies strongly with effective 
temperature and surface gravity, and less strongly with 
metallicity. Drake (1991) finds that the non-LTE effect on 
abundances of Ca in G and K dwarf stars increase con- 
siderably with decreasing effective temperature. From his 
Figs. 4, 7 and 8, we estimate the correction factor for weak 
lines in a dwarf star of at least solar metallicity with an 
effective temperature of 4500K to be on the order of 0.3 
dex. Such an adjustment would indeed put the K dwarf 
stars right on the line, [Ca/Fe] = 0.0 dex. This suggests 
that the calcium abundance may vary in lockstep with the 
iron abundance also for metal- rich K dwarf stars. 

6.8. Titanium 

[Ti/Fe] was shown by Edvardsson et al. (1993a) to be a 
slowly decreasing function of [Fe/H]. The decline may con- 
tinue also for higher iron abundances. We use 10 — 12 
lines to derive titanium abundances for our stars, while 
Edvardsson et al. (1993a) used four. In spite of our, pre- 
sumably, smaller random errors in the abundance determi- 
nation for each star, as is shown in Fig. [ujwe still find the 
same and comparatively large scatter in the abundances 
found by Edvardsson et al. (1993a). 

Inspection of derived stellar abundances as a function 
of excitation energy for the lower level in the transition 
for each line indicated no presence of non-LTE effects or 
blends. However, also for titanium we found no evidence 
that the large star-to-star scatter should be a result of a 
mixing of stars with different mean-perigalactic distances, 
i.e. different Vlsr velocities, see Fig. ^. 

6.9. Scandium and Vanadium 

Abundances derived from Sc I lines are unreliable and we 
therefore only present abundances determined from Sc n 
lines, Fig. [l9| Scandium exhibits some scatter but seems 
to vary in lockstep with iron. 

For vanadium the atom is represented by two lines and 
the ion by one. We present the data derived from lines of 
the atom. Also vanadium appears to vary in lockstep with 
iron over the metallicity range studied. 

6.10. Chromium 

For the most metal-rich stars chromium, as well as other 
iron peak elements, varies in lockstep with iron, Fig. || 
The overall evolution of [Cr/Fe] seems to be well described 
by Timmes et al. (1995) using the supernova yields by 
Woosley & Weaver (1995), Fig. |o| The flatness of the 
relation between [Cr/Fe] and [Fe/H] can be understood 
as a consequence of that massive stars with solar initial 
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Fig. 20. Chromium data from several sources, as given in the figure, showing the galactic chemical evolution of 
chromium. The solid line is from Timmes et al. (1995). The chromium rich star in our sample is HD87646. The error 
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et al. (1993a). 



metallicities produce enough chromium to balance the iron 
production by SNIa, cf. Timmes et al. (1995). 

6.11. Manganese and Cobalt 

Few studies of stellar abundances have been made for 
these elements. Manganese abundances were measured 
by Gratton (1989) for 25 metal-poor giants and dwarfs. 
Cobalt abundances were obtained by Gratton & Sneden 
(1991) for 17 metal-poor (mostly) giant and dwarf stars 
and by Ryan et al. (1991) for 19 dwarf and giant stars. 



Five lines of the manganese atom were used for de- 
termination of abundances, three weak and two stronger 
lines. [Mn/Fe] scales with [Fe/H], but with a tendency 
to increase for [Fe/H] > —1.0 dcx. This increase seems to 
continue beyond [Fc/H] > 0.0 dex. We use seven lines aris- 
ing from the atom to determine cobalt abundances. From 
our data cobalt seems to vary in lock-step with iron for 
0.0 dcx < [Fc/H] < 0.3 dcx. 

Manganese and cobalt belong to the group of iron- 
peak elements. These elements are thought to be formed 
during explosive silicon burning in supernova explosions 
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and nuclear statistical equilibrium (Woosley & Weaver, 
1995). 

Timmes et al. (1995) find that the rise in [Mn/Fe] vs 
[Fe/H] from [Fc/H] ~ — 1 dex is due to the over-production 
of manganese in supernovae type la and in heavy stars 
with solar metallicity, as compared with the iron produc- 
tion, while for cobalt the production of iron in supernovae 
type la is balanced by production of cobalt in supernovae 
resulting from massive stars with initially solar metallicity. 



6.12. Nickel 

We have used 12 lines from the nickel atom to obtain 
abundances; Edvardsson ct al. (1993a) used 20 lines. Like 
Edvardsson et al. (1993a) we find that nickel varies in lock- 
step with iron, and this continues also for higher metal- 
licities. For the stars in common between the studies the 
nickel abundances derived are in excellent agreement, al- 
though our study does not show the slight offset found by 
Edvardsson et al. (1993a). 

Here, we also note an interesting behaviour of the K 
dwarf stars, namely that they show larger nickel abun- 
dances than the rest of the sample. The large number of 
lines together with the fact that the lower excitation ener- 
gies for the lines span a range of values (1.9 — 5.3 eV) and 
that the formal error for each star is small makes depar- 
tures from excitation equilibrium an unlikely explanation 
for this effect. 

Overionization, for most stars but less for the cool 
ones, or blends are possible but neither very probable 
explanations. The phenomenon needs further systematic 
study. 



6.13. r- and s-process elements 



Most of the heavy elements (A>70) are formed through 
the r- and s-processcs. For some of them one of the pro- 
cesses contributes much more than the other. The s- 
proccss contributes most, for the solar system composi- 
tion, to Y (73%) and Zr (79%) while Eu is to 90% formed 
in the r-process, according to Anders & Grevesse (1989). 
Eu is one of the few r-process elements with clean lines 
observable in the visual part of stellar spectra. Therefore, 
it is well suited for studies of the sites for the r-process. 
The relative abundances of s-elements produced in ther- 
mally pulsing asymptotic giant branch stars are set by 
the degree of the exposure to neutrons. Heavier neutron 
flux enhances the abundances of the heavier elements (Ba, 
Nd, Hf) relative to the lighter ones (Y, Zr, Mo). Molyb- 
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Fig. 22. Comparison of data from this study and from 
Edvardsson et al. (1993a) for three elements, yttrium, 
zirconium and neodymium. o symbols denote our abun- 
dances, x symbols denote the K dwarf stars, * symbols 
the stars from Barbuy & Grenon (1990) and • symbols 
denote abundances from Edvardsson et al. (1993a). 



denum is formed by a mixture of processes (p- , r- and 
s-processes), see Anders & Grevesse (1989). 

From this knowledge one would expect the r-process 
elements to have high abundances in old stars and show 
a declining trend when compared to iron. This is indeed 
seen for europium, see figures and discussions in Mathews 
et al. (1992) and Woolf et al. (1995). For the s-process 
elements, on the other hand, one would expect old stars 
to have low s-element abundances while the more recently 
formed stars would show an increase in their s-elemcnt 
abundance, due to the long time scales for the evolution 
of the s-process sites. Such a tendency was also traced by 
Edvardsson et al. (1993a). 
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Fig. 23. Molybdium, lanthanum and hafnium abundances 
relative to iron, x symbols denote K dwarf stars and * 
symbols the stars from Barbuy & Grenon (1990). Note 
the different ranges on the vertical axes. 



Results Our stellar spectra contain few lines of these 
heavy elements that are accessible in the visual and may 
be securely used as abundance criteria. 

We have derived abundances for a number of s- and r- 
process elements, using a small number of Yn lines of suit- 
able strength, one line each for Zn and Moi, two for La 11 
and one each for Ndn, Euli and Hfn. We find that the 
metal-rich stars have roughly solar abundances of these 
elements relative to iron, however, with some possible de- 
parting trends. 

Up to [Fe/H]~0.2 dex we confirm the result by Ed- 
vardsson et al. (1993a) that Y varies in lock-step with iron. 
For the more metal-rich stars, however, there may be a de- 
cline in [Y/Fe] (see Fig. ^]). The results may, however, be 
due to possible effects of overionization in yttrium. 

As is clear from Fig. for zirconium there is probably 
a systematic trend with effective temperature, resulting in 
a large scatter (or even a division of stars into groups) and 
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Fig. 24. The sums of light and heavy s-process elements 
are shown. The sums has been weighted as follows light = 
0.5 ■ [Y/Fe] + 0.3 • [Zr/Fe] + 0.1 • [Nd/Fe], 
heavy = 0.25 ■ [Mo/Fe] + 0.5 ■ [La/Fe] + 0.25 • [Hf/Fe]. 
Only stars with at least two of the elements in the sum 
measured are shown, o symbols denote the sum of the 
light s-process elements and x symbols denote the heavy 
elements. 



unreliable abundances. Blends may may also contribute to 
this scatter. The results for our stars do generally indicate 
lower neodymium abundances than those of Edvardsson et 
al. (1993a) and again a decrease of [Nd/Fe] with increasing 
[Fe/H]. " " 

For molybdenum, lanthanum and hafnium we have 
no studies to compare with and can therefore not say 
very much about the general evolution. The molybdenum 
abundances are derived from one Moi line. Mo and Hf, 
and probably also La, however, show the familiar pattern, 
ascribed to overionization in the K dwarf stars. 

In conclusion, the abundances of Mo, La, and Hf seem 
to roughly vary in lock-step with Fe, however, with some 
indications that the abundance ratios decrease with in- 
creasing [Fe/H]. 

In order to improve the statistics we have derived two 
quantities, < a -light/ Fe >= 0.5- [Y/Fe] + 0.3- [Zr/Fe] + 
0.1 • [Nd/Fe], and < s - heavy /Fe >= 0.25 • [Mo/Fe] + 
0.5 • [La/Fe] + 0.25 • [Hf/Fe]. The weights in these ex- 
pressions reflect the number of spectral lines measured of 
each element. The results are plotted vs [Fe/H] in Fig. 
|24| , where the K dwarfs have been excluded. A downward 
slope of roughly the same magnitude for the "light" and 
the "heavy" elements is seen. This trend is in fact also 
present in most of the corresponding diagrams for the in- 
dividual elements, although the scatter in larger, see Figs. 
|22| and ^3|. A slope of this magnitude can also be traced 
for the s-elements in Edvardsson et al. (1993a), Fig. ^2|. 

A tendency of this type may, if true, indicate that s- 
elements enrichment occurs less frequently in metal-rich 
AGB stars. One may speculate that this might be because 
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Fig. 25. Europium abundances, o symbols this work, x 
symbols K dwarf stars this work, • symbols Woolf et 
al.(1995). Iron abundances are derived from Fei and Fen, 
in the two panels respectively. 

mass loss could finish their evolution earlier than for more 
metal-poor stars. 

6.14- Europium. 

In some studies the iron abundance derived from Fe n lines 
are preferred as reference element for europium rather 
than abundances derived from Fei lines. From our data 
this does not seem to be an obvious choice. Particularly, 
this is not so for the K dwarf stars, since europium with its 
low ionization energy, 5.7 eV, will remain highly ionized 
for all our stars irrespective of the effective temperature 
while almost all iron will be in the neutral state in the 
cool stars. Thus abundances derived from Fe n for these 
stars would be vulnerable to departures from LTE in the 
ionization equilibrium. In our study we derive iron abun- 
dances from, in general, more than 30 lines from Fe I and 
from four or three lines from Fen. Thus, also statistically 
we could expect the atomic abundance to be better deter- 
mined. 

In Fig. ^5] data from Woolf et al. (1995) and our data 
are plotted with iron abundances derived from Fe I and 
Fe II as reference, respectively. Europium shows a declining 
trend with metallicity from [Fe/H] = —1.0 to 0.0 dex and 
this trend now seems to continue unchanged for [Fe/H] 
> 0.0 dex. 

Europium is well-correlated with oxygen as well as 
with the a-elements, Fig. |2(| This supports the idea that 
europium, oxygen and the a-elements are all formed in 
the same type of events, supernovae type II. 



28 



Sofia Feltzing & Bengt Gustafsson: Abundances in metal-rich stars 



0.5 I 1 1 

x ..-<* 

X..-' 
X ° 

o V 
' o 

^ 0.0 - • . ° 

ID 
UJ 



-0.5 - 



-0.5 0.0 0.5 

[0/H] 

Fig. 26. Europium abundances and oxygen abundances 
compared, o this work (excluding the K dwarf stars), x 
K dwarf stars, • Woolf et al.(1995). The line with slope 
+ 1 is indicated by a dotted line. Oxygen abundances are 
in our work derived from the [O i] line at 6300 A while 
the oxygen abundances used together with the europium 
data from Woolf et al.(1995) are from Edvardsson et al. 
(1993a) and are derived (mostly) from the triplet lines at 
7774 A, scaled to the [O i] abundances using results from 
Nissen & Edvardsson (1992). 



It should finally be mentioned that the abundance 
trends (or lack of trends) discussed here are consistent 
with the results obtained for 9 solar-type dwarf stars in 
the interval 0.05 dex < [Fe/H] < 0.25 dex by Tomkin et al. 
(1997). One exception from this, however, is Eu for which 
the latter results sooner suggest a slight increase in the 
[Eu/Fc] with increasing [Fe/H] on the basis of measure- 
ments of the same Eu lines as used here. 



7. Conclusions Summary 

We have determined elemental abundances for 47 G and 
K dwarf stars with iron abundances ranging from solar to 
[Fe/H] =0.4 dex. The stars selected either move on rela- 
tively eccentric orbits spending most of their time inside 
the solar circle, even having been born there, or move 
on solar like orbits. Among the stars on solar like orbits 
there are also a number of Extreme Population I stars, 
with photometric logarithmic metallicities [Me/H] > 0.4 
dex. These samples were selected in an attempt to study 
differences in the chemical evolution in different parts of 
the galactic disk. Our main results and conclusions are 
summarized below. 

The results fall roughly into two categories; galactic 
chemical evolution as studied from G dwarf stars and is- 
sues related to spectroscopy of dwarf stars and in partic- 
ular to ovcrionization in K dwarf stars. 



Galactic Chemical Evolution Our data do, in general, fit 
well into the overall picture of galactic chemical evolution 
as it is currently understood from empirical studies. They 
give further constraints on theoretical models of galactic 
chemical evolution for the region above solar metallicity, 
where metallicity dependent supernova yields further com- 
plicate the picture. In particular we find: 

1. that Mg, Al, Si, Ti, Ni, Ca, and Cr follow iron also at 
these high metallicities, 

2. that [Fe/H] docs not vary with Vlsr, nor any other 
element. This indicates that the chemical evolution 
has lately been very homogeneous in the regions under 
study. This is in contrast with findings for a-elemcnts 
in the more metal-poor part of the disk population, 
Edvardsson et al. (1993a). 

3. that oxygen continues to decline with increasing 
[Fe/H]. This has also been confirmed by Castro et al. 
(1997). 

4. that the upturn in [Na/Fe] vs [Fe/H], as observed 
by Edvardsson et al. (1993a), is reproduced and rein- 
forced. The origin of the upturn is, however, still some- 
what unclear, but metallicity dependent supernova 
yields certainly play a role. In contrast to sodium wc 
find no such upturn for aluminium. This is somewhat 
surprising since sodium and aluminium arc thought 
to be produced in the same environment in the pre- 
supernova star. 

Our lack of success in tracing variations in relative 
abundances, e.g. of [O/Fe], [Na/Fe], [Si/Fe], [Ti/Fe] etc 
to vary with the stellar velocities, i.e. with the present or- 
bital mean distance from the galactic centre (i? m ). Thus, 
the probable mixture of the metal- rich stellar sample, e.g., 
of one old and one more recent population is not revealed 
as differences in relative abundances. That may be because 
of several reasons. One is that the relative number of stars 
from each population may not change very strongly with 
galactic radius. Another may be that the yields of heavy 
elements from the sites of nucleosyntheis are such that, 
e.g., the interstellar Na/Fe ratio may stay rather constant 
through the evolution of the galactic disk. A third possi- 
bility is that galactic orbital diffusion is strong enough to 
bring these differences below the limit of detection. Thus, 
the observed scatter in such abundance ratios can not be 
explained as the result of the mixing of stars with differ- 
ent formation sites at different distances from the galactic 
centre. 

For oxygen and europium our data suggest a contin- 
uation in the declining trend found by earlier studies at 
lower metallicities. This suggests that enrichment from su- 
pcrnovae la is still important, also at these relatively high 
metallicities. 

The scatter in titanium abundances found in Edvards- 
son et al. (1993a) is reproduced in our data in spite of 
the fact that we use more lines to determine abundances 
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from. It seems probable, but remains to be proven, that 
this scatter is "cosmic" . 

Stars with very high metallicities, so called Super 
Metal Rich, SMR, stars, have attracted some attention 
over the years. These stars have been claimed to either 
represent a late stage of the local galactic disk evolu- 
tion or to be members of the Galactic bulge, which have 
now been scattered out to the solar neighborhood. Bar- 
buy & Grenon (1990) obtained photometric metallicities 
and oxygen abundances for eleven such stars. They found 
these stars to have high photometric metallicities and to 
be over-abundant in oxygen. Castro et al. (1997) have ob- 
tained abundances for 9 similar stars. They find that oxy- 
gen, calcium and titanium relative to iron declines with in- 
creasing [Fe/H] , while sodium, silicon and nickel relative to 
iron increases. Several of these stars have fairly low tem- 
peratures and we interpret their results for calcium and 
nickel as further evidence for the phenomenon possibly 
arising from overionization. No velocity data have so far 
been published for these stars and therefore it is difficult 
to compare with our results. We have, however, obtained 
spectra and derived abundances for three of the stars in 
Barbuy & Grenon (1990). These stars remain puzzling - 
their proposed very high [Fe/H] have not been verified, 
but at least some of them seem to have, as these authors 
suggest, a high [O/Fe] ratio. This should be further ex- 
plored. 

Overionization For the K dwarf stars we do detect what 
appears to be a case of overionization for, in particular, 
iron and chromium, such that abundances derived from 
Fe II and Cr II are much higher that those derived from 
Fe I and Cr I. Fe I abundances do not vary with the effective 
temperature. Also for calcium and nickel the K dwarf stars 
depart markedly from the trend outlined by the G dwarf 
stars. Signs of similar departures for rare earth atoms are 
also traced. 
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Table 1. Elemental abundances our programme stars. For each star and each ion the derived elemental abundance, the line- 
to-line scatter (the error in the mean = line-to-line scatter/^/n lines ) and the number of lines used in the analysis are given. The 
table is continued on the following pages. 



ID [Oi] Nai Mgi All Sii Cai Sci 



HD30562 






0.21 


0.06 


(2) 


0.33 0.14 


(3) 


0.25 


0.01 


(2) 


0.21 0.06 


(3) 


0.17 0.09 


(4) 




HD36130 


0.21 


(1) 


0.16 


0.01 


(2) 


0.04 0.24 


(2) 


0.18 


0.03 


(2) 


0.10 0.03 


(3) 


0.11 0.11 


(4) 




HD37088 


0.29 


(1) 


0.16 


0.04 


(2) 


0.21 0.20 


(3) 


0.11 


0.11 


(2) 


0.18 0.07 


(2) 


0.10 0.05 


(4) 




HD37216 






-0.07 


0.03 


(2) 


0.07 0.13 


(3) 


0.01 


0.02 


(2) 


-0.05 0.16 


(2) 


-0.05 0.06 


(4) 


0.19 (1) 


HD49178 


0.16 


(1) 


0.13 


0.00 


(2) 


0.11 0.07 


(3) 


0.00 


0.01 


(2) 


0.00 0.14 


(3) 


-0.02 0.03 


(3) 




HD54322 


0.33 


(1) 


0.18 


0.01 


(2) 


0.11 0.03 


(2) 


0.21 


0.08 


(2) 


0.12 0.04 


(3) 


0.21 0.05 


(4) 


0.48 (1) 


HD55693 






0.34 


0.01 


(2) 


0.43 0.10 


(2) 


0.30 


0.01 


(2) 


0.29 0.12 


(3) 


0.19 0.10 


(3) 




HD67228 


0.21 


(1) 


0.23 


0.01 


(2) 


0.23 0.04 


(3) 


0.22 


0.04 


(2) 


0.27 0.12 


(3) 


0.15 0.02 


(3) 




HD68988 






0.49 


0.10 


(2) 


0.31 0.08 


(2) 


0.37 


0.00 


(2) 


0.39 


(1) 


0.36 0.04 


(3) 


0.69 (1) 


HD69582 


0.14 


(1) 


0.01 


0.01 


(2) 


0.30 0.45 


(2) 


0.07 


0.03 


(2) 


-0.01 0.08 


(3) 


0.02 0.08 


(4) 




HD69830 


0.05 


(1) 


0.01 


0.01 


(2) 


0.00 


(1) 








-0.03 


(1) 


-0.11 0.09 


(3) 


0.26 (1) 


HD71479 


0.24 


(1) 


0.31 


0.02 


(2) 


0.32 0.15 


(3) 


0.29 


0.00 


(2) 


0.30 0.11 


(3) 


0.20 0.08 


(4) 




HD72946B 






0.22 


0.01 


(2) 


0.24 0.27 


(3) 


0.23 


0.02 


(2) 


0.13 0.08 


(3) 


0.18 0.11 


(4) 


0.54 (1) 


HD75782 






0.31 


0.02 


(2) 


0.31 0.26 


(2) 


0.26 


0.07 


(2) 


0.38 0.06 


(2) 


0.22 0.07 


(2) 




HD76780 






0.25 


0.06 


(2) 


0.10 0.34 


(2) 


0.26 


0.01 


(2) 


0.14 0.05 


(3) 


0.13 0.11 


(3) 




HD80607A 


0.41 


(1) 


0.52 




(1) 


0.56 0.15 


(3) 


0.34 


0.04 


(2) 


0.45 


(1) 


0.17 0.09 


(3) 


0.32 (1) 


HD87646 






0.55 


0.05 


(2) 


0.20 0.12 


(2) 


0.35 


0.02 


(2) 


0.28 0.11 


(2) 


0.36 0.04 


(2) 




HD91204 






0.31 


0.01 


(2) 


0.31 0.12 


(2) 


0.22 


0.02 


(2) 


0.31 0.01 


(2) 


0.19 0.04 


(2) 




HD94835 


0.04 


(1) 








0.28 0.06 


(2) 


0.10 


0.03 


(2) 


0.18 0.03 


(2) 


0.10 0.00 


(2) 




HD101242 


0.13 


(1) 


0.05 


0.03 


(2) 


-0.08 0.13 


(2) 








-0.05 0.01 


(2) 


0.03 0.10 


(4) 


0.40 (1) 


HD106156 






0.24 


0.00 


(2) 


0.43 0.10 


(2) 


0.19 


0.03 


(2) 


0.12 0.04 


(3) 


0.21 0.08 


(2) 


0.42 (1) 


HD110010 


-0.04 


(1) 


0.53 


0.01 


(2) 


0.48 0.12 


(3) 


0.40 


0.06 


(2) 


0.51 


(1) 


0.34 0.07 


(3) 




HD 11 7243 


0.13 


(1) 


0.47 


0.06 


(2) 


0.35 0.06 


(2) 


0.32 


0.03 


(2) 


0.39 0.08 


(2) 


0.24 0.14 


(4) 




HD125968 






0.18 


0.01 


(2) 


0.13 


(!) 


0.13 


0.01 


(2) 


0.20 


(1) 


0.13 0.06 


(2) 


0.50 (1) 


HD126511 


0.18 


(1) 


0.16 


0.01 


(2) 


0.16 0.04 
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